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Abstract 

This  paper  presents  results  from  the  self-consistent  solution  of  Schrodinger 
and  Poisson  equations  obtained  in  one-dimensional  non-uniformly  doped 
MOS  structures  suitable  for  the  fabrication  of  very  short  transistors.  Dif¬ 
ferent  issues  are  considered  and  investigated,  including  quantum-induced 
threshold  voltage  shifts,  low-field  electron  effective  mobility  and  gate-to- 
channel  capacitance.  The  reported  results  give  indications  for  the  opti¬ 
mization  of  n-MOS  channel  doping  profiles  suitable  for  the  fabrication  of 
ultra-short  MOSFETs. 


199711)2  054 


Introduction 


As  it  is  well  known,  the  MOSFET  scaling  theories  (e.g.  [1,  2])  give  detailed  prescriptions 
on  how  both  physical  dimensions  and  technological  parameters  should  be  changed  to 
properly  design  short  channel  devices.  In  particular,  high  doping  concentrations  (ss 
1018cm-3)  and  thin  gate  oxides  («  3nm)  are  required  to  limit  2D  short-channel  effects 
(SCE)  in  ultra-short  (US-)  MOSFETs  with  gate  length  (Lg)  below  0.1  /j, m  [3], 

At  the  Si-SiC>2  interface  of  MOS  systems,  electron  energy  quantization  occurs.  In 
the  case  of  the  earliest  MOS  technologies,  its  main  consequence  consisted  in  the  reduc¬ 
tion  of  low-field  mobility  with  respect  to  bulk  Silicon  [4].  In  modern  technologies,  the 
effects  of  quantization  tend  to  become  more  relevant  and  easily  detectable  [5,  6,  7,  8]. 
In  particular,  high  channel  doping  concentration  is  responsible  for  strong  energy  quan¬ 
tization  that,  in  turn,  causes  quantum  induced  threshold  voltage  shift  (Q-DV)  [5,  6]. 
Furthermore,  as  the  centroid  of  the  charge  distribution  is  displaced  from  the  interface, 
gate  capacitance  reduces  with  respect  to  C ox  =  cox/tox  [7,  8],  and  such  a  reduction 
becomes  more  relevant  as  t ox  is  scaled  down  and  the  channel  doping  is  decreased. 

Therefore,  theoretical  studies  based  on  semi-classical  models,  aiming  at  the  opti¬ 
mization  of  the  channel  doping  profile  [9],  must  be  extended  to  account  for  the  effects 
of  energy  quantization.  Previous  works  have  been  devoted  to  the  analysis,  on  a  quan¬ 
tum  basis,  of  the  threshold  voltage  deviations  caused  by  high  dopings  [5,  6],  and  to  the 
effects  of  quantization  on  the  inversion-layer  capacitance  [7,  8].  However,  the  consid¬ 
ered  devices  did  not  feature  extremely  non-uniform  channel  doping  profiles,  like  those 
obtainable  by  epitaxial  techniques  [10]. 

In  this  study,  a  self-consistent  one-dimensional  (ID)  Poisson-Schrodinger  solver  [11] 
is  used  to  determine  the  electrostatic  behavior  of  devices  featuring  highly  non-uniform 
doping  profiles,  which  are  more  realistic  for  the  fabrication  of  US-MOSFETs.  As  a  case 
stud}',  quantum  effects  in  an  epitaxial-channel  (EPI)  MOS  structure  are  investigated. 
Since  we  are  interested  in  the  operation  of  MOSFETs  at  room  temperature,  all  the 
simulations  are  performed  at  300  K. 

To  compare  the  different  channel  doping  profiles,  three  macroscopic  quantities  are 
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evaluated:  i)  the  quantum-induced  threshold  voltage  shift;  ii)  the  electron  low-field 
mobility;  iii)  the  MOS  capacitance,  including  the  effect  of  charge  displacement  from 
the  interface. 

The  obtained  results  put  in  evidence  that  quantum  phenomena  in  MOS  inversion 
layers  are  not  simply  related  to  doping  values  close  the  Si-Si02  interface,  but  a  more 
complex  relationship  holds,  involving  the  actual  doping  profile  within  a  range  of  few 
tenths  of  a  micron  away  from  the  interface.  Furthermore,  this  study  provides  design 
guidelines  for  the  optimization  of  channel  doping  profiles,  to  be  traded-off  with  those 
related  to  2D  SCE. 

The  paper  is  organized  as  follows:  Section  1  presents  the  simulated  device  structure, 
the  adopted  physical  models,  and  the  numerical  techniques  used  in  this  study;  Section  2 
presents  the  simulation  results.  Finally,  Section  3  summarizes  the  conclusions  of  this 
work. 


1  Simulated  structure  and  methodology 

In  this  work,  we  simulated  a  one-dimensional  MOS  capacitor  featuring  an  ideal  dop¬ 
ing  profile,  as  could  be  obtained  by  epitaxial  growth.  A  low-doped  EPI  (with  doping 
concentration  Npp/  and  thickness  tpp/)  covers  a  highly  doped  region  (doping  concen¬ 
tration  Nyj)  acting  as  a  background  plane  to  suppress  punch-through.  The  purpose  of 
the  surface  layer  is  to  improve  low-field  mobility  by  reducing  the  transverse  electric  field 
(TEF)  (i.e.  surface  roughness  scattering)  and  the  effect  of  ionized  impurities.  Since  we 
are  investigating  a  possible  structure  for  MOS  devices  with  Lg  <  0.1  fi m,  we  assumed 
very  thin  oxide  (typically  t0x=3nm)  and  epitaxial-layer  thicknesses  (ts/  as  thin  as  10 
nm).  Fig.l  reports  the  diagram  of  the  conduction  band  across  the  MOS  structure  and 
the  first  20  quantized  energy  levels  self-consistently  computed  for  the  case  toA'=3nm, 
tEP/=10nm  N^=lxl018  cm-3,  Npp/=1  xlO15  cm-3,  inversion  electron  sheet  concen¬ 
tration  Ns=3xl012  cm-2. 

As  said,  the  issues  of  Q-DV,  low-field  electron  effective  mobility  and  gate-to-channel 
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capacitance  have  been  investigated  by  a  one-dimensional  approach.  Therefore,  the 
results  obtained  are  valid  for  low  applied  drain-to-source  voltages,  i.e.  in  the  linear 
MOSFET  regime,  or,  more  in  general,  when  two-dimensional  effects  can  be  assumed 
to  be  small.  Despite  to  this  limitation,  the  results  of  this  paper  are  anywav  valid 
as  a  comparison  between  the  different  doping  solutions  analyzed.  In  fact,  since  the 
concept  of  ’’short  MOSFET”  (i.e.  suffering  of  relevant  2D  short  channel  effects)  is 
relative  to  the  adopted  technology,  even  a  0.1  micron  MOSFET  fabricated  with  a  50 
nm  technology  and  operating  at  low  voltages  [3],  can  be  considered  as  a  ’’long  channel” 
device.  Therefore,  for  such  a  device  ID  results  are  relevant. 

The  simulations  were  carried  out  by  self-consistently  solving  the  Poisson  and  Schrodinger 
equations  with  a  Gummel-like  iteration  [11].  The  envelope  function  equation  (i.e. 
Schrodinger  equation  in  the  effective  mass  approximation)  is  solved  to  determine  the 
eigenvalues  and  eigenvectors  of  the  system.  In  this  solution,  the  six-fold  ellipsoidal 
symmetry  is  assumed  for  the  silicon,  with  the  usual  values  for  longitudinal  and  trans¬ 
verse  masses  (0.19  and  0.915  in  free  electron  mass  units,  respectively),  together  with  a 
parabolic  energy  vs.  k-vector  dispersion  relationship,  adequate  for  the  case  of  n-channel 
MOS  structures  [12].  A  zero  wave  function  boundary  condition  is  forced  at  the  quan¬ 
tum  system  boundaries,  i.e.  in  the  substrate  and  at  the  oxide  interface.  Therefore, 
wave  function  penetration  into  the  oxide  is  neglected. 

From  the  numerical  point  of  view,  the  solution  is  determined  in  the  frame  of  the 
finite  elements  formalism,  solving  the  generalized  eigenproblem  using  the  inverse  iter¬ 
ation,  Rayleigh-quotient,  and  Gram-Schmidt  algorithms.  - 

Once  the  eigenstates  are  determined,  the  quantum  electron  density  and  the  clas¬ 
sical  hole  density  are  computed  using  Fermi-Dirac  and  Maxwell-Boltzmann  statistics, 
respectively.  Complete  ionization  of  the  dopant  is  also  assumed. 


4 


2  Simulation  Results 


In  this  Section  we  report  the  simulation  results  related  to  energy  quantization  inside 
the  channel  of  an  EPI  device.  For  the  sake  of  clarity,  the  discussion  has  been  orga¬ 
nized  into  three  subsections,  dealing  with  threshold  voltage,  effective  mobility,  and  gate 
capacitance  effects,  respectively. 

2.1  Threshold  voltage  results 

We  have  simulated  the  EPI  structure  with  t epi  down  to  10  nm  and  t ox  down  to  3nm, 
investigated  in  [9]  for  application  to  MOSFETs  with  Lq  down  to  50  nm.  Its  threshold 
voltage  (Vth)->  computed  as  the  linear  extrapolation  of  the  simulated  quantum  inver¬ 
sion  charge  Q5  vs.  voltage  characteristic,  has  been  studied.  In  particular,  we  focused 
on  the  threshold  voltage  shift  with  respect  to  the  classical  solution  (i.e.  the  Q-DV) 
for  the  EPI  structure  and  for  the  more  conventional  highly  doped  (UHD)  structure 
obtained  setting  t#p/~0.  Due  to  the  low  surface  doping,  a  much  lower  Q-DV  could 
be  expected  in  the  EPI  case  compared  to  UHD.  Instead,  similar  Q-DV  is  obtained 
for  UHD  and  EPI  with  small  t epi-  In  fact,  for  t^p/— lOnm  and  t ox  ranging  from  3 
to  10  nm,  we  obtained  only  a  20%  reduction  of  Q-DV  with  respect  to  UHD,  showing 
that  this  effect  does  not  simply  depend  on  the  doping  at  the  interface,  but  rather,  a 
long-range  dependence  on  doping  holds.  This  point  is  enphasized  by  Fig.2-a,  where 
Q-DV  is  reported  as  a  function  of  t epi-  a  relatively  thick  EPI  is  needed  to  eliminate 
the  quantum  effects  originated  by  the  presence  of  the  high-doped  background  plane. 
Fig.2-b,  instead,  shows  how  Vth  can  be  controlled  by  properly  selecting  t epi'  while 
the  UHD  MOS  presents  a  threshold  voltage  that  is  too  large  for  low  supply  voltage 
applications  (V dd  ~  1-1.5  V)  the  introduction  of  the  low-doped  layer  provides  a  further 
degree  of  freedom,  additional  to  the  doping  of  the  background-plane,  for  designing  the 
threshold  voltage  according  to  circuit  applications.  The  reported  results  clearly  show 
that  simple  Q-DV  models  based  on  average  channel  doping  values  cannot  be  applied 
to  the  highly  non-uniform  cases  needed  for  US-MOSFETs. 
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2.2  Effective  mobility  results 

We  can  address  channel  electron  effective  mobility  starting  from  Fig. 3.  showing  the 
profiles  of  the  quantum  electron  density  in  structures  that  feature  X.4=lxl0ls  cm-3, 
Nep/=1  xlO  cm  3,  tox=3nm,  and  different  t epi-  The  comparison  between  differ¬ 
ent  structures  has  been  performed  for  a  given  inversion  sheet  density  above  threshold 
(Ns=3x1012  cm  2).  For  increasing  t epi  charge  confinement  reduces  due  to  the  reduc¬ 
tion  of  transverse  electric  field.  Consequently,  an  improved  effective  mobility  (peff  ) 
shall  be  expected  for  larger  t epi- 

In  fact,  the  introduction  of  low-doped  EPI  effectively  decreases  the  effective  electric 
field  (EEF),  reported  in  Fig. 4  as  a  function  of  t epi  for  a  given  above-threshold  inversion 
sheet  density.  Here,  the  EEF  has  been  computed  as  the  mean  value  of  the  electric  field 
component  normal  to  the  interface  averaged  over  the  carrier  densitv. 

The  corresponding  /j,e^  ,  extracted  from  universal  mobility  curves  [4|  is  also  given 
in  Fig.4,  confirming  what  anticipated.  It  should  be  mentioned  that  the  UHD  device 
is  expected  to  present  even  larger  disadvantages  for  low  inversion  charge  densities  [4], 
due  to  the  reduced  screening  of  the  Coulomb  scatterers. 

Additional  results  (Fig.  5)  show  that  a  lower  effective  field  is  obtained  for  the 
EPI  with  respect  to  the  UHD,  also  when  the  comparison  is  performed  for  the  same 
gate  drive  (e.g.  Vg-Vt'h),  which  is  a  very  significant  comparison  in  view  of  circuit 
applications. 

2.3  Gate  capacitance  results 

Another  important  issue  relates  to  total  AfOS  capacitance  Gtot-i  i-e.  the  senes  of  the 
oxide  capacitance  Cox=£ox/tox  and  of  the  inversion  layer  capacitance  C  ixv=dQs/dtps, 
tps  being  the  Si-surface  potential.  As  mentioned,  C tot  in  quantized  inversion  layers 
reduces  with  respect  to  Gox  as  t ox  is  decreased,  due  to  the  larger  effective  oxide 
thickness  consequent  to  the  displacement  of  the  inversion  charge  from  the  device  inter¬ 
face.  As  shown  in  [7]  for  uniform  channels,  such  a  reduction  becomes  more  relevant  for 
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lower  dopings.  Consequently,  compared  to  the  UHD  case  a  lower  capacitance  should 
be  expected  for  the  EPI  one.  Fig.6-a  reports  Ctot—^-Qs / ^ g  as  a  function  of  t ox-,  for 
^EPI~  10  nm  and  Onm.  As  can  be  seen,  the  capacitance  degradation  increases  at  thin¬ 
ner  oxides,  but  no  additional  degradation  occurs  if  thin  low-doped  EPI  is  introduced. 
As  reported  in  Fig.6-b,  the  simulated  results  indicate  that  no  serious  degradation  shall 
be  expected  by  the  application  to  US-MOSFETs,  such  as  those  described  in  [9],  of  EPI 
layers  up  to  30-40  nm.  Results  obtained  by  classical  simulation  accounting  for  Fermi- 
Dirac  statistics  are  also  reported.  Degenerate  statistics  can  only  partially  account  for 
the  capacitance  attenuation  occurring  at  thin  gate  oxides  as  already  reported  in  [8], 
furthermore,  classical  simulation  underestimates  capacitance  dependence  on  t epi- 

3  Conclusions 

In  this  work  a  self-consistent  Schrodinger-Poisson  simulator  has  been  applied  to  the 
study  of  highly  non-uniform  doping  profiles  suitable  for  the  fabrication  of  ultra-short 
MOSFETs.  The  obtained  results  prove  that  a  non  trivial  dependence  of  quantum  ef¬ 
fects  on  the  channel  doping  profile  holds.  The  introduction  of  a  low  doped  region  at 
the  device  surface  allows  to  tailor  the  threshold  voltage  according  to  circuit  applica¬ 
tions,  where  the  use  of  high  dopings  would  lead,  instead,  to  threshold  values  too  high 
for  realistic  applications.  Furthermore,  the  presence  of  a  low-doped  epitaxial  region 
produces  a  reduction  of  the  electron  effective  field  for  a  given  charge  sheet  density 
and,  therefore,  improves  the  electron  effective  mobility.  As  for  gate  total  capacitance, 
the  displacement  of  the  inversion  charge  from  the  silicon  surface  causes  its  reduction 
with  respect  to  Cox,  but  no  significant  additional  degradation  has  been  found  as  a 
consequence  of  the  introduction  of  a  low-doped  epitaxial  layer. 
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Figure  Captions 

Fig-1  Conduction  band  across  the  MOS  structure  Fermi  level  (dash-dotted  line)  and 
the  first  20  quantized  energy  levels  computed  for  the  case  to*  =3  nm,  tEPl~10  nm 
Ny4=lxl018  cm-3,  N epi  =  lxlO15  cm-3,  inversion  electron  sheet  concentration 
Ns=3xl012  cm-2.  The  Fermi  energy  is  only  7  meV  below  the  first  eigenvalue. 

Fig. 2  a)  Quantum  induced  Vjh  shift  vs.  EPI  thickness;  b)  Vth  vs.  EPI  thickness. 
to*=3nm,  N#p/=1015cm-3,  N,i=1018cmr3. 

Fig. 3  Electron  density  vs.  position  for  different  EPI  thicknesses.  Symbols  are  for  vi¬ 
sualization  only  and  do  not  correspond  to  actual  discretization  points.  The  bias 
Vg  was  set  in  order  to  obtain  the  same  charge  sheet  density  for  all  structures. 

Fig-4  a)  Effective  electric  field  (squares,  left)  and  effective  mobility  (circles,  right)  vs.  EPI 
thickness.  The  bias  VG  was  set  in  order  to  obtain  the  same  charge  sheet  density 
for  all  structures. 

Fig-5  Effective  electric  field  for  the  uniform  MOSFET  (circles)  and  for  the  EPI  MOS 
with  tg;p/=30  nm  vs.  Gate  Drive  (Vg-Vx//). 

Fig- 6  a)  Total  capacitance  vs.  oxide  thickness  for  tpp/=0,  10  nm.  Symbols:  simulated 
results.  Dashed-line:  Co*,  b)  Total  capacitance  vs.  t epi  for  to*=3  nm.  Circles: 
simulated  results;  crosses:  results  of  classical  simulation.  Dashed-line:  Co*. 

All  simulated  structures  feature  Npp/=1015cm~3,  N^=1018cm“3.  The  bias  VG 
was  always  set  in  order  to  obtain  the  same  charge  sheet  density  (Ns=3x  1012cm-2) 
for  all  structures. 
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